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B. F. Robinson and M. E. Bauer 
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Laboratory for Applications of Remote Sensing 
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ABSTRACT 


In this paper we present the results of an experiment de- 
signed to characterize a soybean field by its reflectance at 
various view and illumination angles and by its physical and 
agronomic attributes. Reflectances were calculated from mea- 
surements at four wavelength bands through eight view azimuth 
and seven view zenith directions for various solar zenith and 
azimuth angles during portions of three days. An ancillary 
data set, consisting of the agronomic and physical character- 
istics of the soybean field is described. These data sets 
should prove useful for validating most light interaction 
canopy models. 

The results of the study indicate that the distribution 
of reflectance from a soybean field is a function of the solar 
illumination and viewing geometry, wavelength and row direc- 
tion as well as the state of development of the canopy. 

Shadows between rows greatly affected the reflectarce in the 
virible wavelength bands and to a lesser extent in the near 
infrared wavelengths. 

A model is proposed that describes the reflectance vari- 
ation as a function of projected solar and projected viewing 
angles. The model appears to approximate the reflectance 
variations in the visible wavelength bands from a canopy with 
well defined row structure. 


1. INTRODUCTION 

The reflectance characteristics of an agricultural scene are a function of 
agronomic, geometric and atmospheric variables. Accurate characterization of 
agricultural crops from remotely measured spectral radiance or spectral reflec- 
tance is necessary if these data are to be used effectively for identification 
and condition assessment purposes. 

In the past, researchers have noted that identification of agronomic vari- 
ables from remote sensing data is affected by the dependence of canopy reflec- 
tance on the solar illumina' ion and sensor viewing geometry. It is well known 
that increasing our understanding of this phenomenon will improve our abilities 
to monitor the physiological and phenological status of a crop. 


•Presented at the Fifteenth International Symposium on Remote Sensing of Envi- 
ronment, Ann Arbor, MI, May 1981. 
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In recent years two goals have pron¥>ted researchers to measure the reflec- 
tance properties of crop canopies for various illumination and view conditions. 
First, the knowledge gained from analysis of en^iricai data, primarily acquired 
for nadir viewing conditions, has increased our understanding of spectral re- 
flectance properties as they relate to the physiological and phenological sta- 
tus of a crop canopy. Daughtry et al , 1980, found relatively high correlations 
between spectral reflectance and such agronomic variables as biomass, percent 
ground cover, leaf area index and development stage of wheat. Learner et al, 
1980, showed that reflectance could be used to monitor crop development. A 
strong relationship between spectral radiance and grain yield of wheat was des- 
cribed by Tucket et al, 1980. Holben et al, 1980, found significant correla- 
tions between spectral radiance and soybean leaf area and biomass. All or the 
above studies utilized either truck-mounted or hand-held nadir viewing sensors. 
The need to extend this understanding to off-nadir viewing conditions is in^ior- 
tant due to the variations in view direction of airborne and satellite multi- 
spectral scanners and the selective pointing capabilities of future satellites 
such as the Systeme Probatoire d’ Observation de la Terre (SPOT) (Chevrel et al, 
1980) to be launched in this decade. 

The second goal for acquiring reflectance data at off-nadir viewing angles 
is to provide data suitable for testing the accuracy of mathematical models 
which have been proposed to describe the radiance of a plant canopy for speci- 
fic view and illumination conditions. An accurate model of the plant canopy 
reflectance could potentially provide the same information as is obtained from 
field observations, but without extensive field data acquisition activities. 
Analysis results of reflectance measurements have demonstrated that plant cano- 
pies arc not simple Lambertian reflectors and thus not amenable to description 
by simple models (e.g. , Suits, 1972; Smith and Oliver, 1974; Rimes et al, 1980; 
Vanderbilt, 1980). Since the models are not simple, their veracity must be 
tested vigorously, using data sets acquired under a variety of view/illumina- 
tion and crop physiological and phenological conditions. As most canopy re- 
flectance data have been acquired for nadir viewing conditions, only limited 
data acquired under a wide variety of view/illumination conditions are available 
for testing models. No extensive data set exists for soybeans, a crop of major 
economic importance. 

In this paper we present the results of an experiment designed to charac- 
terize a soybean field by its reflectance at various viev; and illumination an- 
gles and by its physical and agronomic attributes. We discuss the reflectance 
differences obtained over a variety of view and illumination angles as affected 
by differences in the shadow between rows caused by the relationship between 
illumination, viewing and row directions. In addition, we propose a simple mod- 
el to explain those reflectance variations. This model incorporates the con- 
cepts of projected solar angle and projected view angle to explain a portion of 
the reflectance variation due to shadows cast by rows. 

2. MATERIALS AND METHODS 

A commercial soybean field located approximately 16 km northeast of West 
Lafayette, Indiana was selected for this study. The field was planted to Cala- 
tian 9250 soybeans at a seeding rate of 62 kg/ha in north-south oriented rows 
planted 76 cm apart, resulting in about 28 plants/n/. The field is in corn- 
soybean rotation and was hand-weeded prior to data acquisition to remove any 
volunteer corn and other weeds. 

Reflectance Measurements 


The reflectance measurements were acquired on three clear days in the sum- 
mer of 1980 (July 17 , July 24 and August 28) with an Exotech model 100 radio- 
meter with wavelenuth bands of . 5- . 6 , .6-. 7, ,7-, 8 and .8-1.1 ’,m. The field of 
view of the instrument was limited to lO' by a set of field stops attached to 
the viewinu ports. Tho instrument was positioned on a truck-mounted aerial 
platform at a nominal elevation of 10 meters above the soil surface. The in- 
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struRient was mounted on a panhead which alisjwed it to be positioned at specific 
view zenith and azimuth angles. The truck and associated instrument van and 
electrical generator were driven into the middle of the soybean field suf£i> 
ciently distant from field edges to uniformly fill t)w field of view of the 
radiometer. 

The measurement procedure consisted of positioning the boom of the aerial 
platform at an azimuth angle (from north) of 135” (Figure 1). At this position 
the instrument was set at a view zenith angle (6v) With held con- 

stant, spectral radiance measurements were made for view azimuth angles (8v) of 
0" (looking south), 315°, 270® (looking east), and 225”. This procedure was 
repeated for Oy = 15”, 22”, 30”, 45” and 60”. The boom was rotated 180” and 
the process was repeated for = 180”, 135”, 90” and 45* with By “ 7”, 15”, 

22”, 30”, 45” and 60”, thus completing a hemisphere of measurements. In addi- 
tion to the off-nadir measuren^ts, vertical canopy radiance was acquired at 
boom azimuth angles of 90®, 135”, 160”, 270”, 315” and 340”. In order to convert 
spectral canopy radiance to spectral reflectance factor (Robinson and Biehl, 
1977), a 1.3 X 1.3 m calibration panel was viewed before and after each hemi- 
sphere measurement sequence. Just prior to the last calibration measurement, 
an area cleared of soybeans was measured with nadir view angle to provide data 
for bare soil spectral reflectance. Color photographs were taken with a boom- 
mounted 35mro can«ra coincident with each scene measurement. These photos were 
later used to identifv any measurements where a boom shadow might influence the 
calculated reflectance. 

This procedure was repeated at roughly one-half hour Intervals as long as 
cloud conditions remained satisfactory (Table 1) . Care was taken to avoid tak- 
ing moasurements when the sun was obscured or when individual cumulus type 
clouds were within 15" of the solar disk. 

Meteorological data consisting of relative humidity, air temperature, 
barometric pressure, wind direction, wind speed and global solar irradiance 
were measured on each day at the Purdue Agronomy Farm located about 6 km south- 
east of the study site. As an indication of the irradiance variation at the 
site, an upward- looking photocell was monitored throughout the n^asurement 
periods. In addition, shaded calibration panel measurements were taken on each 
of the three days to estimate the proportion of skylight present. 

Agronomic Measurements 

Agronomic and physical characteristics of the soybean field were measured 
the day after each reflectance data acquisition date. The agronomic measure- 
ments consisted of sampling fresh and dry biomass, leaf area index, average 
canopy height and width and development stage. Physical measurements included 
leaf inclination angle, leaf azimuth angle, leaf area and three-dimensional 
leaf position. In addition, individual leaf reflectance and transmission niea- 
surements were obtained using live plants. 

Soybean biomass was estimated by harvesting plants in 1.0 m lengths of row 
at randomly located positions within the field. Each sample was placed in plas- 
tic bags and weighed to determine fresh biomass. The soybeans were then sepa- 
rated into green leaflets, yellow leaves, stems plus petioles and pods, dried 
at 6'’” C and reweighed. A random subsample of green leaf area from each sample 
was measured with an optical device (Lambda Instrument Co. Model 1-3000) and 
the ratiii of leaf area to leaf dry weight was estimated. Leaf area index was 
calculated using this leaf area: Leaf weight ratio and the total dry weight of 

green leaflets from plans in a 1.0 m length of row. Development stage was es- 
timated based on the Fehr and Caviness scale (Fehr and Caviness, 1977). In ad- 
dition, estimates of percent ground cover were made from the vertical color 
photographs obrained during reflectance data acquisition. 
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Physical Measureiwnts 

’mrae'dlmansional characterisation of the soybean canopy was aceoaplished 
by measuring leaflet heights, distances fr<»n mid-row and distances along row as 
references to a 1 m x 1.3 m board with a grid of xaeter sticks attached to it. 
Leaflet inclination angles were determined with an inclinometer similar to those 
described by Nichlporov^ch (1961) and Kyle and Davies (1974). Asimuths were 
measured with a simple apparatus consisting of a snail plate equipped with a 
circular bubble level and a dowel attached perpendicular to the face of the 
plate. 

The location of the sanple was randomly determined by throwing a meter 
stick out into the field. The measurement board was then placed in this loca- 
tion in the furrow halfway between adjacent rows and a windbreak was constructed 
to protect the sanple from wind. The board was then leveled and positioned ver- 
tically. The distance between the bottim of the board and the base of the 
plants was measured to correctly determine the leaflet heights. All leaflets 
within a 25 cm length of row were sanpled starting with the uppermost leaf and 
progressing downward. Leaflet azimuth angles were measured by orienting the 
bubble level plate towards true north and then placing it parallel to the upper 
leaflet surface. The bubble level was graduated into 45° segments with the po- 
sition of the bubble indicating the azimuth + 180°. With the azimuth device 
still in place, the inclinometer was held parallel to the perpendicular dowel 
on the plate, thus indicating the angle from vertical to the leaf normal. The 
recorded angle was later converted to inclination. After the angle measure- 
ments were made, a meter stick was placed above the leaf and the distance from 
the center of the leaflet to the measurement board was recorded. Being sure 
the meter stick was perpendicular with the measurement board, we recorded the 
distance from a reference point on the board to the center of the leaflet and 
the height of the leaflet. The leaflet was then removed from the plant, num- 
bered with a felt marking pen and placed in an ice chest for later leaf measure- 
ment. Individual leaf areas were determined with an optical leaf area meter 
later the same day in the laboratory. 

The geometry estimation procedure was performed simultaneously by two 
teams at different locations in the field. Each crew, consisting of a measure- 
ment person and a recorder, received over 4 hours training prior to the first 
data collection period. All destructive sampling was conducted outside the 
total field of view perimenter for the reflectance measurements. 

3. RESULTS 

The estimates of the agronomic variables for the three days are summarized 
in Table 2. The soybean canopy developed rapidly between the first and second 
measurement dates with ground cover increasing from 72% to 83%, indicating that 
the canopy was increasing in density and rapidly losing its previously well de- 
fined row structure. By the last measurement date, the canopy was completely 
overlapping with a ground cover of 99%; however, LAI had decreased about 10%. 

The distribution of leaf inclination angles changed rapidly from a more or 
less planophile type (mostly horizontal leaves) on the first date to having a 
greater proportion of vertically inclined leaves on the second and third dates. 
Blad and Baker (1972) reported that soybean leaves distribution tends to be 
more planophile through most of the growing season, with leaves distributed al- 
most evenly in the early and late growth stages. We suspect that the departure 
of the leaf distribution on the second day from a planophile type to a more 
evenly distributed arrangement is due to wilting caused by the very high temper- 
atures and lack of rainfall prior to and during this date. Figure 2 presents 
the leaf angle distributions weighted by leaf area for the three dates. 

To evaluate how the changes discussed above affect the reflectance from 
the soybean field, reflectance factors for each of the four wavelength bands 
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were plotted over eolar zenith angle, 6 ^ £> 3 r a nadir sensor viewing angle (dy “ 
0°) (Figure 5). Ilte soybean canopy reflectance factor (Rq) increased slightly 
with increasing for the two near infrared wavelength bands. Quite a differ* 
ent trend was observed for the ttio visible bands for the first two dates. 
was the largest in both bands when the sun was at its highest position in the 
sky, then decreased to a point and then leveled off through steeper 6 ^. This 
trend was reported by Kollenkark et al, 1981, to be mainly the effect of (hang- 
ing shadows cast between rows of soybeans. The point at which the reflectance 
begins leveling off is a function of 6 g,^g and the canopy ^ometry. fffhere ca- 
nopy geometry is described in terms of the across-vow canopy profile (Figure 3), 
it can be shown that a critical solar angle ( 6 ^) defines the point where soil 
between the rows seen by a nadir looking sensor becomes completely shaded. Thus 
at 63 < sunlit bare soil is present in the scene. 

The reflectance factory in the visible band for the third measurement date 
when the canopy was overlapping follows the same trend as the near IR reflec- 
tances. This suggests that completely developed canopies are affected to a les- 
ser extent oy changes in 63 than canopies with well defined row structure. This 
isn't meant to imply that overlapping soybeans approximate a Leunbertian surface; 
rather, quite the opposite is true, as we will discuss later. 

The reflectance of bare soil decreased gradually as 65 increased for all 
four wavelength bands for all three dates (Figure 5). This gradual decrease 
was not due to shadows cast by soybean rows, but was more likely due co surface 
roughness effects caused by soil clods and furrows. 

It is apparent from the above results that the reflectance in the two visi- 
ble wavelength bands of a soybean field is strongly affected by shadows cast by 
rows. This effect is not as strong in the near infrared bands nor is it appar- 
ent in the visible bands for an overlapping canopy. 

Analysis of the reflectance factors for varying view and illumination geo- 
metry shows quite dramatically the non-Lambertian character of the soybean field. 
Figure 6 presents polar graphs of equal reflectance contours for the . 6 -. 7 ym 
and . 8 - 1.1 pm wavelength bands over a hemisphere of viewing zenith and azimuth 
angles at two different times for each date. To interpret these graphs, the 
reader is referred to Figure 4 which describes the positions of the view zenith 
and azimuth ancles of the measurement hemispheres. 

For the July 17 measurements, graphs are presented for hemispheres of re- 
flectance factors acquired at midday and late afternoon (Table 1) . In the first 
hemisphere, the distribution of reflectance factors for the visible band can be 
likened to a hill elongated along the row direction ($-). Reflectance falls 
off more rapidly at i 90° than it does at ^r 1 180°. The near in- 

frared reflectance, however, is maximum in the direction toward the sun and de- 
creases to a minimum away from the sun at Oy = 30°, and then increases slight- 
ly in the same direction for steeper Oy. For the second hemisphere where 63 is 
at nearly right angles with 4>j., the reflectance factors approximate a bowl with 
a steep side towards the sun. A similar pattern occurs with the near infrared 
band reflectances. 

For July 24, hemispheres are presented for midmorning and midday. The vi- 
sible reflectance distribution for the first hemisphere shows a somewhat simi- 
lar pattern to that of the first July 17 hemisphere with a maximum occurring 
hear 6 y = 20°, ty = 0°. The near infrared reflectance is highest at Oy * 30° 
in the direction towards the sun and decreases away from the sun. The second 
hemisphere has maximum visible band reflectance at 9y » 40° in the direction of 
the sun and a minimum at 0 „ = 20° away from the sun. The near infrared reflec- 
tance is also maximized towards the sun at Oy « 50° and decreases away from the 
sun. The August 27 hemispheres acquired at midday and late afternoon all show 
bowl -like distribution with the steepest side toward the direction of the sun 
for both wavelength bands. 
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It appears that the reflectance distribution of the soybean field has a 
maximum where the viewing and illumination angles coincide near the hot spot. 

We can infer this only indirectly since when spectral data were acquired at the 
hot spot, the field of view contained the shadow of the instrument and these 
data were not included in this analysis. 

In the .6-. 7 ym wavelength, region shadows between rows are the dominant 
feature in the scene. The amount of shadows is a function of the canopy geo- 
metry and the position of the sun relative to the row. Reflectance measurements 
from a vertical- looking sensor are a function of the changing shadows as the 
solar zenith angle changes. As the view angle changes, more or less shaded 
soil and canopy are viewed, thus either decreasing or increasing the reflectance. 


If the sun azimuth is in line with the row direction, then the amount of 
shaded soil will be minimized. With the view azimuth rise equal to the row di- 
rection, the canopy reflectance will be maximized at the hot spot. If the view 
azimuth direction is at right angles to the sun and row direction, reflectance 
will decrease as the view zenith angle increases due to less sunlit soil and 
nwre canopy shadows being present. 

If the sun's azimuth is at right angles tr row direction, then reflectance 
will be maximized at the hot spot and will decrease at smaller view angles due 
CO shaded soil. If the view azimuth is 180“ from the sun, reflectance increases 
with view angle, possibly due to forward-scattered light. 


we can describe the relationsnip of viewing and illumination angles with 
row direction as projected view angle (9vp) and projected solar angle (6sp)» 
discussed by Verhoef and Bunnik, 1976, and Koilenkark et al, 1981. That is, 

'* = TAN~l(tan6s sin (Is “ 'I’r) ) and ^’vp = TAN~^(tan6y sin(ly - • Addition- 

ally, we can describe a critical angle {Qq) as that Sgp w.iere the soil is 
completely shaded. 


Let us first assume a square-shaped canopy where only sunlit foliage and 
sunlit soil contribute to the canopy radiance (i.e., shaded foliage and soil 
are black) and the sunlit portion of the canopy approximates a Lambertian sur- 
face. Then we can describe a model that should account for varying reflectance 
as a function of view and illumination geometry. We identify limits of Ogp, 

'vp and 9^ where the reflectance should have nearly equal reflectance. When 
^0'’ 'so c ^VD '"'sp' have the case of the canopy hot 

spotT fn this region only raSianbe from the sunlit top of the canopy in^inges 
the sensor and R(0gp, ‘Vp) will approximate R(9gp, ) . For the case where 


90 


_ sp _ 


and 


0 


— - r — - vp .. - - 

ow so R(0gp, tiyp) will be less than the hot spot reflectance and equal to the 


sp » 


the sensor sees proportionately more shad- 


case where 0‘ 


sp 1 "c 0 ' ^vp ' ^^sp' Assigning a negative sign to view 

directions from the left side of the row^and making the additional assumption 
that the canopy is symmetrical along the row, then when 90° 2. ^sp ' -90* 

■ '\,p ■ 0“or Hg •< 'V and -90“ • t'Vp ' 0“, then R(9gp, O^p) will be about 


'sp _ c “““ vp “ ' viio ■'v'spf '^vp' 

to the reflectance for a given and a nadir-looking sensor. If we 


and . iVp 1 9sp' then if the canopy is viewed from 


cqua 

have _ 

the reflectance will be least when 9vp = *V and will increase up to the hot 
spot position. Finally, when 0“^ -sp ■ “c and 0“ ■ 'Vd ' "sp' the sensor 

will see the same amount or shaded soil and canopy ana the reflectance will ap- 
proximate R(9gp, ’’c) . A summary of these results is presented in Table 3. 


Areas of equal reflectance cap be identified on the basis of the differ- 
ence between Ogp and «vp' as shown in Figure 7. Analysis of our reflectance 
data in the . 6- . 7 ;,m band indicates this model is a reasonable first approxima- 
tion to the complex distributions of reflectance obtained over a wide range of 
view and illumination angles. It is, however, quite apparent that it fails to 
deal with the multiple scattering, both diffuse and specular, that complicates 
the reflectance distribution in both the visible and near-infrared bands. 
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4. StMNJUty 

The rather conq>lex distribution of reflectance from a soybean field over 
varying view/ illunina ton conditions has bMn discussed. It is evident that the 
reflectance distribution is a function of the solar and viewing ge«Mtry« wave- 
length, and row direction as well as the stare of i.«velopMnt of the canopy. 
Shadows between rows greatly affected the reflects e fron the field in the 
visible wavelength bands and to a much lesser exten'.. the near-infrared band 
reflectances. 

The results of the experiments discussed should provide a data set suitable 
for testing various canopy models that predict the reflectaiuie as a fui^tion of 
view and illumination conditions, canopy gec^etry and agronomic variables. A 
document further describing this data set will be forthcoming later this year. 

A model was proposed that describes the expected variation of reflectance 
in terms of the difference between projected solar and projected view angles. 

The assumptions of the model limits its applicability to row canopies and the 
highly absorptive red wavelength band. It does, however, provide an initial 
fran^worlc for understanding the rather coRf>lex reflectance variations of a soy- 
bean canopy as a function of view and illu^nation conditions. 
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Table 1. Zlluisinatior. Conditions for ^ctral Data Collactim 


solar 




Solar zenith 



Angle 


Tine Period 

Angle hettee 

Solar Aiinuth 

ClOttd 

Period 


Start 

Stop 

Hex. > Nin. 

• 

Max. 

Anole Range 

Cover 





Degrees 


Degrees 

1 

Datet July 17> Day of Yean 199 






1* 


17tS9 

18:32 

19 

- 

21 

183-20S 

10 

2 


19t03 

19t29 

24 

- 

28 

222-233 

10 

3 


19t47 

20t06 

31 

- 

34 

240-245 

20 

4 


20i28 

20t49 

38 

- 

42 

251-256 

10 

5* 


21tl9 

21:35 

47 


50 

262-265 

10 

Datet 

July 

IS, Day 

of Year: 206 





1* 


15tl4 

15:31 

40 - 37 



109-113 

1 

2 


IStSl 

16:06 

33 - 31 



119-124 

20 

3 


16t31 

16:51 

27 - 25 



133-142 

10 

4 


17tl0 

17:25 

23 - 22 



152-161 

5 

5* 


17tSl 

18:C6 

21 



178-188 

5 

6 


18t36 

18:49 

22 

- 

24 

206-214 

20 

Datet 

August 27, Day of Yeart 

240 





1 


15:55 

16:14 

40 - 37 



132-139 

0 

2 


16t23 

16:39 

36 - 34 



142-148 

0 

3 


17t07 

17:23 

32 - 31 



160-167 

0 

4 


17: 30 

17:45 

31 - 30 



171-178 

0 

5* 


18tl7 

18:31 

31 

- 

32 

193-200 

0 

6 


18t3S 

18:52 

32 

- 

33 

203-209 

0 

7 


19t06 

19:20 

35 

- 

36 

214-220 

0 

8 


19:27 

19:43 

37 

- 

39 

222-227 

0 

9 


20:10 

20:23 

43 

- 

46 

235-239 

0 

10 


20:32 

20:46 

47 

- 

49 

241-245 

0 

11* 


21:08 

21:26 

53 

- 

56 

249-253 

0 

12 


21:42 

21:59 

59 


60 

256-237 

0 

* Indicates 

measurement hemispheres shown in Figure 6. 







Table 2. 

Agronomic Measurements 



Date 


Leaf 

Area 

Index 

Total Fresh 
Biomass 
g/nr 

Total Dry 
Biomass 
q/n* 

Ground 
Cover % 

Height 

cm 

Width 

cm 

Maturity 

Stage 

July 

18 

3.0 

1145 

2 30 

72 

69 

58 

V13R3 

July 

25 

3.9 

1540 

320 

83 

84 

69 

V14R3 

Aug. 

28 

2.8 

2535 

645 

99 

102 

104 

V20R6 
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Figure S. JRelationships of soybean canopy and bare soil reflectance factors and solar zenith 

angle on three days for a nadir looking sensor. 
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Figure 7. Ntxlel diagreni of e,p plane. 

Numbered regions within solid lines indicate 
areas of nearly constant reflectance defined 
by 6ap< @vp and Oc (see Table 3). Dashed 
lines represent equal values of 
angles. ^ 


Table 3. Sumnary of Reflectance Relationships 

for Angular Limits of © ,6 and 6 . 

sp vp c 

Area of ~ ~ " 

©sp't'vp Angle Limits Reflectance 

Plane 


1 

90* ' r 
- »P 

- ®vp 

-^.p 

«<^p'V 

- «<®sp' V 

a constant 

2 

90® ' e 
- »P 

- «c' 0 ' ®vp ^ 

^■p 

«<®SP'V 

'«‘^P'V 


i 

90° > e 

- sp 

' 0 , -90° • 0 
- c ' vp 

» 0° 

«<^p'V 

a constant a 

R<e,p.0) 

4 

0°- 1* 

- sp 

' i 

®»P 

R(0,p,e^) < 

«‘%p'«vp» 

' "<®sp'®.p 

5 

0° ^ © 

- sp 

«c' % i 

»c 

R(0 ,e„) 

•p vp 

a R(©,p.e^,) . 

« constant 

6 

0® V 0 

- sp 

- ®c' ' ®vp * 

e 

•p 

R(@ .6 ) 

sp' vp' 

^ «‘®sp'V 


7 

^,p 

i °o' ^ ®vp 

< 0° 

R(e ,e ) 
sp vp 


constant 
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